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We calculate the second, fourth and sixth order quark number fluctuations in the deconfined 
phase of 2+1 flavor QCD using lattices with temporal extent N t =4,6,8 and 12. We consider 
light, strange and charm quarks. We use p4 action for valence quarks and gauge configurations 
generated with p4 action with physical value of the strange quark mass and light quark mass 
m q = 0.1m 4 generated by the RBC-Bielefeld collaboration. We observe that for all quark masses 
the quark number fluctuations rapidly get close to the corresponding ideal gas limits. We compare 
our results to predictions of a quasi-particle model and resummed high temperature perturbative 
calculations. We also investigate correlations among different flavor channels. 
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1. Introduction 



Quark number fluctuations are basic observables which can be obtained by taking derivatives 
of the grand canonical potential with respect to quark chemical potentials jj, q . Here we look directly 
at the related derivatives c\ , which enter into the Taylor expansion of the pressure at finite chemical 
potentials (jl q /T < 1) 
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Their combinations form fluctuations of conserved charges (see ref. [[[]] for a related lattice 
study). Thus they can be used to study the properties of a thermal medium produced in heavy 
ion collisions From Eq. ( |1.1[ ) it is clear that the expansion coefficients are related to the 
quark number densities rif, quark number susceptibilities %{ = ( n 
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2<?2 and other moments: 



In this work we are interested in the behavior of the quark number fluctuations at high temper- 
atures 240MeV < T < 900MeV and zero chemical potential. We would like to study the approach 
of relevant observables to the ideal gas (Stefan-Boltzmann (SB) limit) ^ 
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We use the tree level Symanzik improved gauge action and improved staggered fermion p4- 
action with 3-link smearing (p4fat3). The gauge configurations were generated for calculation 
of the equation of state with (2+1) flavors using almost physical light quark masses resulting in 
m n rs 220MeV and physical s quark mass fl]. In our measurements we add non-dynamical c quark. 
In table [l| we list the gauge configurations used in the study. 
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Table 1: Lattice parameters and the number of trajectories for each parameter set. The step size between 
measurements is 100 trajectories for stared entries (*) and 50 for the rest. Typically we start measurements 
after 1000 trajectories. Rnd is the number of random vectors used in the stochastic estimators. 
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2. Numerical results 

Here we present our numerical results for quark number fluctuations up to the 6th order, in- 
cluding second order flavor off-diagonal fluctuations. In Fig. [T] we plot the off-diagonal coefficients 
c"f and c"j normalized by the related diagonal quadratic terms c\ and c\. The plot indicates that 
the correlations between u and d, and u and s flavors vanish rapidly after the deconfinement. This 
is suggestive of the fact that there are no bound states in the deconfinement phase [g]]. Similar 
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Figure 1: The off-diagonal coefficient c\ "j (left) and c\\ (right) normalized by c\ and rfj respectively. 

results were obtained for uc correlations. Thus we may expect that quark gas description should be 
reasonable for quark number fluctuations. 
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Figure 2: Numerical results for c\ calculated for different N t and the corresponding SB values (left). The 
ratio c s 2 /c l 2 as function of the temperature (right). 



In the left panel of Fig. ^ we present our lattice results for c\ coefficients as well as the 
corresponding continuum SB limit. In addition we plot the SB value corrected for a finite lattice 
cutoff. In Fig. |2] we also show our previous (low temperature) results from Ref. [|l|]. We see 
that the observed cutoff effects are different and smaller than the cutoff effects in the free theory. 
The difference between the N t = 6 and N t = 8 results is about 10%. Therefore to establish the 
continuum limit calculations on Af f = 12 lattices are needed. Similar results were obtained for c s 2 . 
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The right panel of Fig. || shows the ratio c\jc\. As this has been discussed in |[TJ, |6|] this ratio 
approaches unity at T > 300MeV. This in fact is expected in the quark gas picture. 

In the left panel of Fig. || we compare our c\ results for T > 200MeV with the resummed high 
temperature perturbative calculations This figure clearly shows the 10% difference between 
N t = 6 and N t = & results. Therefore we also did calculations at N t = 12 lattices at 3 temperatures. 
The corresponding results are also shown in the figure. At the highest temperature, where A 7 , = 12 
data are available, we see a reasonably good agreement with N t = & data. This suggests that N t = 8 
calculations are sufficiently close to the continuum results. Therefore the comparison of N t = 8 data 
with the resummed perturbative calculations performed in the continuum limit is meaningful and 
shows that quark number fluctuations are well described by the resummed perturbative calculations. 
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Figure 3: Left panel: The ratio of c\ to the corresponding SB value compared to the high temperature 
perturbative results shown as dashed lines. The three lines correspond to the renormalization scale fi = 
nT,2nT and 4nT (from bottom to top). Right panel: The effective thermal quark mass extracted from c\. 



We also calculated the 4th order expansion coefficients c\ and c\ in the high temperature 
regime. The numerical results for c\ are shown in Fig. || (the results for c\ are similar at high 
temperature). In the figure we show the corresponding SB value in the continuum as well as for 
finite N t . The cutoff dependence of C4 is similar to that of c%. Namely, the observed cutoff effects 
are different from the free theory. Also the ordering of N t = 6 and N t = & results is reversed. The 
difference between them is less than 10%. Assuming that N t = 8 results are close to the continuum 
limit, as the C2 calculations suggest, we see deviations from the ideal gas limit around 20%, i.e. 
two times more than for ci. 

In the resummation scheme used in [j7|] the dominant contribution to the quark number fluc- 
tuations is given by quasi-particles with effective quark mass equal to the hard thermal loop mass. 
Therefore it makes sense to compare the lattice results to a simple massive quasi-particle model. 
The pressure for massive quark gas is given by equation ( |Q| ) H 

£ = f (7) E(- 1 ) m '" 2 ^('-"/n«>*(W), (2.D 
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Figure 4: The high temperature behavior of c\. 
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Figure 5: Numerical results for c% (left) and the effective mass extracted from it (right). 



Taking successive derivatives of the rhs with respect to fl/T we obtain expressions for C2 (|2^2|), 
C4, etc. The expression for C2 can be solved for m/T using C2 values obtained numerically. The 
corresponding quasi-particle masses are shown in Fig. || (right) for N t = 8. We can use these 
effective masses to make predictions for C4 and c&. The quasi-particle results for c ll A is shown in 
Fig. ^ As we can see they are close to the massless ideal gas limit and almost 20% off from the 
numerical data for A 7 , = 8. This indicates that the simplest quasi-particle description is not adequate 
for light quarks. Unfortunately the are no perturbative calculations for c\. It would be interesting to 
see whether the resummed perturbative calculations similar to those performed for c\ in [0] could 
explain the observed deviation from the ideal gas limit for c\. 

Now we consider the charm quark case. In Fig. |5| on the left panel we plot the results for c|. 
Here one can notice that N t = 4 data shows slightly different scaling due to cutoff effects. On the 
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right panel we present the effective mass extracted from the c c 2 values at each temperature point. 
In Fig. ^ we compare c c 4 data with the quasi-particle model results, where the effective mass was 
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Figure 6: Numerical results for c\ and the quasi-particle results. 



computed from c 2 data. We see a very good agreement between the quasi-particle model results 
and the numerical data. This indicates that for heavy charm quark the quasi-particle description is 
adequate for modeling of c c A at high temperature. 

As we see from Figs. |5] and ^ the fluctuations of the charm quark number is comparable 
to those of the light quark sector at the highest temperatures. This means that at temperatures 
T ~ 800MeV the contribution of charm quark to QCD thermodynamics is significant as this has 
been suggested in Ref. ^ using weak coupling analysis. 

We also considered the 6th order coefficient. In Fig. [7] we present the results for c\ and c c 6 . The 



light quark result according to Eq. ( L2j ) should be zero for a massless ideal gas. The numerical 
result agrees reasonably with this prediction. For the charm quark number fluctuation Cg we took 
the effective mass computed from c\ and used it to compute the massive ideal gas estimate, which 
we plot in the figure with lines for various N t . It is interesting that the numerical results for all N t 
are in reasonable agreement with the quasi-particle model results in the deconfined regime. 



3. Conclusions 



In this contribution we discussed quark number fluctuations C2-,c\ and C(, calculated with p4 
action on A 7 , =4,6,8 and 12 lattices. We see that the cutoff dependence of quark number fluctu- 
ations is different from the expectation based on the free theory. Off-diagonal fluctuations vanish 
soon after the deconfinement transition. This suggests that fluctuations can be understood in terms 
of the quark gas. In the light quark sector the quadratic fluctuations are quite close to the massless 
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Figure 7: The sixth order C(, results for c and u flavors. 

ideal gas limit and agree with the resummed perturbative calculations and expectations of the quasi- 
particle model. The quartic fluctuations on the other hand show larger (about 20%) deviation from 
the massless ideal gas limit and do not agree with the predictions of the quasi-particle model. In the 
charm quark sector we find that all fluctuations up to 6th order can be reasonably well described by 
the quasi-particle model. Current results suggest that continuum limit may be reached for N t = 8 
but calculations on ,/V, = 12 lattices at higher temperature are needed to verify this statement. 
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